A model for the production of regular fluorescent light from coherently driven atoms 
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It has been shown in recent years that incoherent pumping through multiple atomic levels provides 
a mechanism for the production of highly anti-bunched light, and that as the number of incoherent 
steps is increased the light becomes increasingly regular. We show that in a resonance fluorescence 
situation, a multi-level atom may be multiply coherently driven so that the fluorescent light is highly 
anti-bunched. We show that as the number of coherently driven levels is increased, the spontaneous 
emissions may be made increasingly more regular. We present a systematic method for designing 
the level structure and driving required to produce highly anti-bunched light in this manner for an 
arbitrary even number of levels. 
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I. INTRODUCTION 



It is well known that resonance fluorescence from a two 
level atom is anti-bunched jl|,||. That is, consecutive flu- 
orescent photons are less likely to arrive very close to- 
gether (ie. to be bunched) than they are in a classical 
light beam of the same average intensity. The distribu- 
tion of waiting-times between photo-detections for reso- 
nance fluorescence is therefore peaked further away from 
zero than for classical light. The reason for this effect is 
easily understood. When the atom spontaneously emits 
to produce a fluorescent photon it places itself in the 
ground state. It cannot emit again until the laser field 
driving the atom has increased the occupation probabil- 
ity of the excited state. The evolution of the excited 
state occupation probability amplitude for a driven two 
level atom (excluding spontaneous emission) is given in 
Fig. (|l|) . It is the nature of this evolution (along with the 
spontaneous emission rate of the excited state) which de- 
termines the form of the waiting-time distribution. 

We note that in recent years it has been shown that 
in laser models with more than two levels, it is pos- 
sible to achieve sub-Poissonian statistics in the out- 
put H . In general sub-Poissonian statistics implies anti- 
bunching p|. Ritsch et al. Q have shown that highly 
sub-Poissonian light may be produced if the atomic state 
is recycled from the lower to the upper lasing level via a 
large number of incoherent pumping steps through inter- 
mediate levels. As the number of pumping steps tends 
to infinity the intensity fluctuations of the output light 
tends to zero. In this paper we will consider a single 
atom resonance fluorescence situation rather than a laser 
model. In this case the driving is completely coherent, 
and the anti-bunched light is produced by a single inco- 
herent transition. 

It has been shown previously that a coherently driven 
three level Raman laser gives a sub-Poissonian output Q . 
Also, Schernthanner and Ritsch |7j have shown that de- 
tuning of the pump and laser light from the atomic tran- 
sitions in a Raman (3-level) laser can lead to a value of 
the Mandel Q parameter close to —1, and hence a highly 



sub-Poissonian output. Ralph and Savage || have shown 
that a coherently driven 4-level laser will produce ampli- 
tude squeezed light. In this paper we treat atoms with 
an arbitrary, but even, number of levels. 

We might call the limit of extreme anti-bunching reg- 
ular light. That is, light in which the time between con- 
secutive photon-detections is constant, or at least, the 
fluctuations in the temporal separation of consecutive 
photons is small compared to the average value of this 
separation. Clearly, if we can arrange the atomic evo- 
lution so that the occupation probability of an unstable 
state makes a sharp transition from zero to some nonzero 
value a well defined time after a spontaneous emission, 
then if we have a high spontaneous emission rate the re- 
sult will be essentially regular fluorescence. 

In the following section we show how we may design a 
coherently driven n level atomic system (where n is even) 
to produce an evolution in which the anti-bunching is in- 
creased over that for a two level system. The evolution 
may be made increasingly closer to a given desired evo- 
lution (for example that to produce re gula r light) as the 
number of levels is increased. In Section III we give an ex- 
ample of this procedure for the case of a four-level system. 
In Section |V] we calculate the waiting-time distributions 
achieved for various n. We conclude in Section |y| 



II. TAILORING THE ATOMIC EVOLUTION 

Excluding spontaneous emission, the equation of mo- 
tion for the occupation probability amplitudes for a 
coherently driven n level atomic system may be writ- 
ten |,0 
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where a,j is the probability amplitude for state j, and H 
is an n x n Hcrmitian matrix, whose elements, we 
will take to be real for simplicity. We are also using units 
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FIG. 1. Solid line: The evolution of the excited state prob- 
ability amplitude of a driven two level atom with Rabi fre- 
quency £1. Dashed line: An imaginary evolution that would 
produce light, the regularity of which would be limited only 
by the spontaneous decay rate of the excited state. 



such that % = 1. The matrix H, which is the Hamilto- 
nian, is determined both by the laser fields driving the 
atomic levels, and the dipole matrix elements between 
the atomic levels. The magnitude of the off-diagonal 
element hij is given by the product of the strength of 
the laser field coupling the atomic levels i and j with 
the magnitude of the dipole matrix element connecting 
those levels. The phase of the off-diagonal element is 
determined by the the phase of the respective coupling- 
laser field and the sign of the respective dipole element. 
The diagonal elements are determined by the detuning 
of the laser fields from the various atomic transitions be- 
ing driven |lf| , [l2]| . The general evolution (up until any 
spontaneous decay) is given by 
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where the coefficients c m are determined by the initial 
conditions chosen for the a,j. The v m are the eigenvec- 
tors of H. We will write their components as v m = 
(vx m , . . . , v nm ), and assume them to be normalised so 
that v m • v ; = 6 m i. 

We now impose certain conditions on the laser interac- 
tion Hamiltonian H . We will show below how interaction 
matrices satisfying these conditions may be constructed. 
Once the interaction matrix has been constructed, this 
tells us the configuration of laser fields that is required to 
drive the atomic levels to produce the desired evolution. 
First of all let us assume that n is even, and that for 
every eigenvalue of H there is another eigenvalue which 
has the same magnitude but opposite sign, so that the 
eigenvalues come in pairs which sum to zero. Let us take 
the initial condition to be a n = 1, and all the other am- 
plitudes zero, so that the atom starts initially in the state 
n. Finally, let us assume that the n values c m v\ m also 
come in opposite pairs, such that c m v\ m — ciVu = iff 
A m — A; = 0. The evolution of a\ may now be written 



which is clearly a sine series for the evolution of a\ . 

Before we show how to construct a matrix H with the 
above characteristics, let us turn to the question of reg- 
ular photon emission. Let us assume that all the levels 
in our system have much lower spontaneous decay rates 
than level 1, and that this level decays via spontaneous 
emission to level n. For regular photon emission we there- 
fore require that the evolution of the amplitude of level 1, 
given by the (finite) sine series above, have a sharp tran- 
sition from zero to some non-zero value. An ideal evolu- 
tion would be given, for example, by the dotted 'square' 
curve in Fig.(|lJ). To approximate the ideal curve to a 
given level of accuracy we may choose the sine series in 
Eq.(||) to be the first n terms in the Fourier sine series 
for the desired 'square' evolution. The Fourier series in 
question is 

* — ' (Zm — 1 7T 

771=1 V ' 

where 1 /ft is the period of f(t). Clearly the evolution will 
more closely follow the desired evolution as n is increased. 

We should note here that we only require the evolution 
to be equal to the Fourier series up to a constant factor. 
This is because the overall scaling of the evolution is not 
important as far as the resulting waiting time distribu- 
tions are concerned; we can cancel the effect of scaling 
the evolution curve by scaling the spontaneous emission 
rate of the unstable level with respect to f2. This will 
be obvious when we write down the expression for the 



waiting time distribution in Section IV 



To construct the required matrix H we first write it in 
the form 



H = TDT 



(5) 



where D is a diagonal matrix containing the eigenvalues 
of H, and T is a matrix whose columns are the normalised 
eigenvectors of H . Constructing H involves choosing the 
elements of D and T. We are clearly free to choose the 
eigenvalues of H, by choosing the diagonal elements of 
D. The second condition imposed above concerns the n 
values c m vi m - Now the c m depend upon the initial con- 
ditions. Recall that we choose the initial conditions to 
be <zi = . . . = a„_i = 0, a n = 1- With this choice it is 
easily shown that c m = v nm . That is, the coefficient c m 
is given by the last element of vector v m ; the condition 
on the coefficients is a condition on the elements of the 
eigenvectors of H . In particular it is a condition on the 
first and last elements of the eigenvectors. 

It is pertinent to note now that if the columns of a 
square matrix are orthonormal vectors, it follows that 
the rows are also orthonormal. To construct H we there- 
fore require both to find a set of n orthonormal rows for 
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T, and satisfy the condition on the first and last elements 
of the vectors v m . The first elements of the v m are the n 
elements of row one, and the last elements of the v m are 
the n elements of row n. The condition we must satisfy 
is that the n values v nm vi m must sum to zero in pairs, 
which clearly implies 



The resulting Hamiltonian is given by 
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This simply states that the dot product of row one with 
row n should vanish. The condition on the elements of 
the v m is therefore consistent with the orthogonality of 
the rows. 

The v\ m and v nrn may now be chosen so that the coef- 
ficients in the sine series for the evolution of a\ match (up 
to an overall constant factor) the desired first n Fourier 
coefficients. Choosing the other n — 2 rows of T is then 
a straight forward procedure in linear algebra. (For ex- 
ample the Gramm-Schmidt orthogonalisation procedure 
may be used [^3[). 

Note that there is clearly a degree of freedom in con- 
structing the matrix H. If we take the matrix T to 
be real, then there are n 2 undetermined coefficients. 
Orthonormality of the rows imposes n(n + l)/2 condi- 
tions. Obtaining the desired Fourier coefficients imposes 
(n — 1) conditions (one less than n because the overall 
scaling of the coefficients is unimportant). This leaves 
(n/2— l)(n— 1) elements undetermined. These may there- 
fore in general be chosen so as to simplify the form of H 
in order to simplify the corresponding physical system. 



III. AN EXAMPLE FOR A FOUR-LEVEL 
SYSTEM 

We now use the procedure described in the previous 
section. We construct a Hamiltonian for a four-level sys- 
tem that will produce evolution corresponding to the first 
two terms in the sine series given by Eq.(^j). 

First we choose the eigenvalues — 67rO, — 27rf2, 27r£l, 
67rfi, which give us respectively the diagonal elements of 
D. To satisfy the condition on the elements of T we first 
choose the top row of T to have elements all identical. 
As a result the bottom row of T must be chosen to have 
elements which vanish in pairs, the absolute value of each 
pair being proportional to the absolute value of a Fourier 
expansion coefficient. Thus we choose as the bottom row 
(1, — 3, 3, l)/%/20- The other two rows are free to be any 
two further mutually orthonormal vectors. We decide to 
choose them such that three off diagonal elements of H 
are zero. The resulting T matrix is 
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H = TDT 
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Note that all the diagonal terms are zero. While the 
off-diagonal elements describe the coupling of the atomic 
energy levels via classical driving fields, the diagonal el- 
ements are determined by the detuning of the driving 
fields from the atomic transitions which they drive. The 
diagonal elements of the Hamiltonian are zero when all of 
the detunings are zero, so that the Hamiltonian we have 
constructed describes the situation where all the lasers 
are tuned to the transitions they drive. The value of 
n determines the strength of the laser fields, and may 
clearly be thought of as a generalised Rabi frequency. 
In Fig. (^j) we give a diagrammatic representation of the 
four-level atomic system showing the coupling between 
the levels. 
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FIG. 2. A diagrammatic representation of the driven 
four-level atomic system designed to produce light with 
anti-bunching increased over a two level system. The ele- 
ments of the Hamiltonian coupling the various levels are given 
as multiples of uj — 27rfi2CT 
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IV. WAITING-TIME DISTRIBUTIONS 

We now derive the waiting time distribution from the 
evolution of the unstable state. Appealing to the Monte- 
Carlo Wave Function (or 'quantum trajectory') method 
for simulating the evolution of open quantum systems, 
the probability for a spontaneous photo-emission in time 

a* is 63| 



7 |ai(t)| 2 At, 



(9) 



where 7 is the spontaneous emission rate of level 1. First 
let P(t) be the probability that there were no emissions 
in the interval [0, t], given that there was an emission at 
time t = 0. The probability that there is no emission in 
the interval [0, t + At] is therefore 
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P(t + At) = P(t)(l- 1 \a 1 (t)\ 2 At). 



(10) 



Waiting-Time Distributions for Various n 



Taking the limit as At tends to zero we arrive at a 
differential equation for P(t), given by 



dP(t) 
dt 



- 7 M*)| 2 p(t), 



the solution of which is 

P(t) = e - 7 i>^')l 2dt ' 



(11) 



(12) 



The waiting-time distribution, w(t), is then P(t) multi- 
plied by the rate of emission at time t. It is clear from 
Eq. (|Tll) that this may be written 



w(t) 



^) =7 | ai ( t0 | 2 e-^ol-(*')l 2d *'. (13) 



dt 



Clearly this equation will provide an analytic solution for 
w(t) as long as the integral of ai(t) possesses an analytic 
solution. Naturally this is true for any truncated Fourier 
series. In Fig.(||) we plot the waiting-time distributions 
that would be generated by systems which produce an 
evolution matching the first n terms in the Fourier series 
given by Eq.([|), for various values of n. We calculate w(t) 
by numerically evaluating Eq.(13), where a(t) is given by 
Eq.(|j). Clearly the waiting-time distribution is increas- 
ingly peaked for increasing n. Experimentally there is a 
great deal of freedom in choosing the generalised Rabi 
frequency f2 since it is determined by the intensity of the 
laser driving fields. For the plots in Fig.(||) we have used 
51 = 7/IOO which is well within experimental limits. 

Another quantity of interest is the second order cor- 
relation function of the fluorescent light, G 2 (t). This is 
the probability density for one photon to be emitted at 
time t and another photon to be emitted at time t + r, 
regardless of how many other photons where emitted in 
the intervening time interval (t,t + r). The second or- 
der correlation function is given by rQ(r), where Q(t) is 
the probability that a photon is emitted at a time t + r, 
given that we have a photon emitted at time t, and r is 
the probability per unit time for a photo-emission. The 
latter is simply the average rate of photo-emission, which 
is 



w{t)t dt 



(14) 



To calculate Q we may use w(t), but we must sum over 
the probabilities for different numbers of photons to be 
emitted in the interval (t,t + r) Jl^| . This is performed 
very neatly by the formula 



Q(r) - w(t) 



Q(t)w(r - 1) dt. 



(15) 




FIG. 3. Waiting-time distributions that would be gener- 
ated by systems which produce an evolution matching the 
first n terms in the Fourier series given by Eq.(^J), for var- 
ious values of n. Solid line: n = 2. Dashed line: n = 4. 
Dash-dot line: n = 8. Dotted line: n = 16. The value of the 
spontaneous emission rate is 7 = lOOfi. 



If we define Q(z) and w(z) as the Laplace transforms of 
Q(t) and w(t) respectively, then this equation may be 
expressed as 



Q(z) = 



w{z) 
1 — w(z) 



(16) 



For a given w(t) we may therefore obtain the correspond- 
ing Q(t), and therefore G 2 (t), by solving Eq.(fL5|) using 
Laplace transforms. Alternatively the second order cor- 
relation function may be obtained by solving the mas- 
ter equation describing the atomic system, and using the 
quantum regression theorem M. For our purposes it is 
simplest to obtain Q(t), and hence G 2 (r), by numerically 
integrating Eq.(|l5|). 

In Fig.(|j) we plot the normalised second order corre- 
lation function, g 2 (r) = G 2 (T)/r 2 , for various values of 
n. As we expect, this has a series of peaks and troughs 
which gradually even out with time. As the regularity of 
the photon emission increases, longer separation times, 
r, are required before the correlation function becomes 
smoothed out. 

We have considered in detail here only the ideal case 
in which spontaneous emission from atomic states other 
than level one can be ignored. However, the effects of 
spontaneous emission from the other atomic levels is eas- 
ily estimated. In this scheme the evolution of the am- 
plitude of level one, ai(t), and the ratio 7/fi are cho- 
sen so that, excluding spontaneous emission from any 
other atomic level, level one will emit after a time de- 
lay of approximately l/(8f2) (that is, approximately one 
eighth through the generalised Rabi cycle which has pe- 
riod 1/51). Clearly we may safely over estimate the ampli- 
tude of the other levels by taking them to be unity. The 
effect of another state with a decay rate of 7 = (80)/fc, 
where k is some dimensionless number, is therefore to 
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FIG. 4. The normalised second order correlation function 
that would be generated by systems which produce an evo- 
lution matching the first n terms in the Fourier series given 
by Eq.(Q), for various values of n. Solid line: n — 2. Dashed 
line: n = 4. Dash-dot line: n = 8. Dotted line: n = 16. The 
value of the spontaneous emission rate is 7 = 100f2. 

interrupt the normal (regular) spontaneous emission se- 
quence on average no more than once every k emissions. 
Taking 7 = 100f2, which is the value we have used for the 
previous plots, then for a value of 7 Q = 7/IO 4 , sponta- 
neous emissions from this level will, on average, interrupt 
the desired sequence fewer than approximately once ev- 
ery 1000 emissions. Clearly in this case, and for larger ra- 
tios k, the waiting time distributions we have calculated 
will be a good approximation to the true distributions; 
any changes to the curves will be on the order of 1/1000 
or less of the height of the main peak. 

Clearly to implement our scheme experimentaly re- 
quires that, for a given atom, atomic levels can be found 
with very different Spontaneous emission rates. Hovever, 
this is not dificult, as spontaneous emmision rates for 
atomic levels are found to vary by many orders of magni- 



tude. Observation of atomic shelving jl6| (interruptions 
of the fluorescence from one 'bright' level by shelving to 
a meta-stable level) with single trapped ions has already 
been performed JL7|. In these experiments the sponta- 
neous emission rate of the meta-stable level is of the or- 
der of 10~ 8 of that of the bright level. An experimental 
realisation of our scheme, at least for the four level sys- 
tem which we have treated explicitly, therefore appears 
to be well within the limits of current technology. 

V. CONCLUSION 

We have shown that multiply coherently driven multi- 
level atomic systems may be designed to produce light 
which is more anti-bunched than for a two level system. 
We have given a procedure for deriving the driving re- 
quired to produce highly anti-bunched light for a given 
number of atomic levels. This procedure requires that 
one of the atomic levels have a spontaneous emission rate 
which is much larger than the spontaneous emission rates 
of the other levels. We have also shown that as the num- 
ber of levels is increased, the fluorescent light may be 
made increasingly regular. 

We have treated explicitly the case of a four-level sys- 
tem, calculating a possible driving configuration to pro- 
duce highly anti-bunched light. 
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